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Abstract — We consider a Gaussian multiple-access channel 
(MAC) with an amplify-and-forward (AF) relay, where all nodes 
except the receiver have multiple antennas and the direct links 
from the transmitters to the receiver can be neglected. Thus, 
spatial processing can be applied both at the transmitters and 
at the relay, which is subject to optimization for achieving 
high data rates. In general, this optimization problem is non- 
convex and hard to solve. While in prior work on this problem, 
it is usually assumed that all transmitters access the channel 
jointly, we propose a solution where each transmitter accesses 
the channel exclusively, using a time-division multiple-access 
(TDMA) scheme. It is shown that this scheme not only simplifies 
the problem a lot, but also yields higher achievable sum-rates. 

I. Introduction 

In today's wireless communication systems, the demand for 
higher data rates and wide-range coverage is steadily growing. 
To meet these requirements, a high density of base stations 
is necessary, which entails high costs for installation and 
maintenance. Another possibility to increase throughput and 
coverage is the use of relay nodes, which have much lower 
costs. Relay channels were considered in [1] first, and have 
drawn more and more research attention in the last decades. 

Depending on how the signals are processed at the relay, 
different types of relaying schemes are distinguished. The most 
common ones are amplify-and-forward (AF, also called non- 
regenerative relaying) and decode-and-forward (DF, also called 
regenerative relaying). While in AF, the relay simply amplifies 
the received signals subject to a power constraint, a complete 
decoding and re-encoding of the signal is necessary when 
using DF. As this yields higher costs and larger delays, we 
will restrict ourselves to AF relaying schemes in this paper. 

In addition to relays, the deployment of nodes with mul- 
tiple antennas helps to further increase the achievable data 
rates. The optimal combination of these two paradigms is an 
enormous challenge and has been considered in numerous 
publications, such as Q-H and the references therein. In 
J2] an optimal amplifying matrix for a MIMO relay is found 
for the case where the transmit covariance matrices of the 
transmitting node is a scaled identity matrix. The general 
case with arbitrary transmit covariance matrices is considered 
in 0, where the optimal structure of both relaying and 
transmit covariance matrix is found. However, this structure 
still contains parameters that are subject to optimization and 
an optimal solution of this problem remains open. 



If multi-user systems are considered, finding the optimal 
transmit covariance and relaying matrices becomes even more 
complex. Achievable sum-rates for broadcast- ||6| 13 and 
multiple-access channels (MAC) [7] with relays have been 
derived, but the optimal rates have not been found yet. 

In this work, we consider a K-user MAC with an AF-relay, 
where all nodes except the receiver have multiple antennas. 
For this system, we first find the optimal solution of the 
approach made in J7). Subsequently, we will introduce a 
different transmission scheme based on time-division multiple- 
access (TDMA), that decomposes the if -user multiple-access 
relay channel (MARC) in K orthogonal single-transmitter 
relay channels. Thus, the relay only receives the signal of 
one transmit station at a time, while, in the transmission 
scheme from [7], the relay receives all transmit signals jointly. 
Comparing the two schemes, we will show that the TDMA- 
based transmit scheme achieves sum-rates that are always 
larger or equal than those achieved by joint relaying of the 
signals in the considered channel. 

This work is structured as follows: In Section [III we intro- 
duce the channel model and describe the constraints that have 
to be fulfilled while optimizing the sum-rate. Subsequently, 
we review a previous solution of this problem in Section [HI] 
and optimize it for the model considered here. A different 
approach, based on TDMA, is proposed in section [IV] After 
optimizing the parameters of the TDMA-based solution, its 
superiority is shown. Subsequently, Section [V] provides a 
quantitative comparison between the solutions by means of 
simulation results. Section [VTl concludes the paper. 

II. Channel model 

A. Notation 

We denote all column vectors in bold lower case and 
matrices in bold upper case letters. The trace and the hermitian 
of a matrix A are identified by tr(A) and A H , respectively. 
We use ||x|| to denote the Euclidean norm of a vector x and 
I to describe the identity matrix. Furthermore, A max (A) and 
Vmax(A) indicate the largest Eigenvalue of a matrix A and 
its corresponding Eigenvector. A diagonal matrix with entries 
oi, . . . , a n is denoted by diag (ai, . . . , a n ). 

B. Channel Model 

The MARC that we consider is depicted in Figure Q] and 
consists of K transmitting nodes. User k G {!,..., K} has 
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Fig. 1. K-usev multiple-access relay channel (MARC) 



M^ k > antennas, transmits the signal x( fc ) E 



relay, and has a channel matrix of g i^M r xM [ - > f or ^ s 
transmission. Thus, the received signal y r of the relay can be 
written as 
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where z r ~ C7V(0, 1) is the additive white Gaussian noise 
at the relay and M r denotes the number of antennas at the 
relay. The relay amplifies the signals by the matrix F and 
transmits the signal x r = Fy r over the channel h H e C lxMr 
to the receiver. It is assumed that the transmission from 
the relay is orthogonal (in time and/or frequency) to the 
transmissions received by the relay and that the direct paths 
between transmitters and receiver can be neglected. Hence, the 
received signal can be written as 

y = h ff x r + z = h H Fy r + z 

= h H F^H«x«+h H Fz r + z, 

fc=i 

where z ~ CJV(0, 1) denotes the Gaussian noise at the 
receiver. Both the transmitters and the relay are subject to 
average power constraints. We use Q( fc ) = E x'*'x''' ) to 



denote the transmit covariance matrix of transmitter k. With 
this definition the transmit power constraints are given by 



tr(Q (fc) ) < P (k) Vfc G {!,..., K} 



E(tr (x r xf )) =tr F I + £ H«Q( fc >H(' 
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k=l 



Finally, we assume that perfect channel state information is 
available at all nodes. 

III. Joint Relaying Scheme 

In Q a transmission scheme was presented, where all 
transmitters send simultaneously and the relay amplifies and 
forwards the sum of all those signals. Therefore, we will refer 
to this scheme as "joint relaying" in the remainder of this 
paper. This scheme is briefly introduced in subsection IIII-AI 
In the following subsections, we first find the optimal relaying 



matrix F (IIII-Bb and then derive the optimal covariance 
matrices Q^, . . . , QW (ITTFCl l. 

A. Previous Results 

The achievable sum-rate of the scheme from [7] is given as 

h H FRF H h 



ife = log 2 1 



h^FF^h+ 1 



(2) 



where R = Ijj^QWH^ . It can be optimized 

with respect to F, Q' 1 ', . . . , Q^. Considering H as a ma- 
trix of size M r x N, where N is the number of receive 
antennas Q, its singular value decomposition (SVD) can be 
defined by H = U^S/jV^ where = diag(<7i, . . . , ctn) 
with descending diagonal order. Moreover, we define the 
Eigenvalue decomposition of R as R = U r A r U^ where 
A r = diag(Ai, . . . , Am,.) with descending diagonal order. 

As stated in [7] and shown in [3], the optimum structure of 
F is given by 

F = U h E / Uf, 
where S / = diag(/ 1; . . . , / M J 1/2 and f % > Vi. 

B. Optimal Relaying Matrix 

For any Q, {1 \ . . . , Q, {K \ we can reformulate (fl} and (O 
with the expressions from above to obtain the following 
optimization problem (cf. Q): 



max i?£ 



min{N,M r } 

E 

i=l 



log 2 



1 



°lk + 1 



s.t. 5^(A< + l)fi < P r and U > Vi 

i=l 

However, for the setup considered here, we have only N = 
1, i.e., <j\ = ||h|| 2 and cr 2 = ■ • • = &n = 0. As only f% 
contributes to the sum-rate, the problem is obviously solved 
by fi = fr(Ai + l) -1 , h = ■ ■ ■ = !mt = for this case. 
Hence, the sum-rate can be expressed as 

cr 2 AiP r 



R s - log 2 ( 1 
= log 2 



a{P r + Ai + 1 



1 



HhfA^R)^ 
h|| 2 F r + A max (R) + 1 

log 2 (l + ||h|| 2 P r ) -log 2 (l 



(3) 



|h|l 2 R, 



1 + A max (R) 

and can now be optimized with respect to Q*- 1 ', . . . , Q^ K \ 

C. Optimal Covariance Matrices 

As it can be seen from the last line of (0, the optimal sum- 
rate is achieved by maximizing the largest Eigenvalue of R. 
As the maximum Eigenvalue operation is a convex function, 
the maximization is not trivial. To provide a solution, we will 
introduce two Lemmas: 

1 H is referred to as h in case of only one receive antenna (N = 1) 



Lemma 1: Let A,B 6 C" x " be positive semidefinite 
matrices. Then, if A — B is positive semidefinite (A >z B), 
we have A max (A) > A max (B) 

Proof: The proof is given in Appendix A. ■ 

Lemma 2: Let II € C nx ™ be a positive semidefinite matrix 
with tr(n) < P and A e C mx ™ . Then, it holds that 

AA H P y AUA H 



channels. Therefore, we will first derive the optimal structure 
of the transmit covariance- and relaying matrix for the sin- 
gle user relay channel in subsection IIV-AI In the following 
subsection IIV-B1 we will transfer this scheme to the MARC 
with TDMA and optimize the achievable sum-rate. Finally, 
in subsection IIV-CI we will show that this sum-rate is always 
larger or equal than the one in ([§). 



Proof: The proof is given in Appendix B. ■ 
With these two lemmas, we can now formulate: 

Theorem 3: The optimal value of the non-convex optimiza- 
tion problem 

max A max (R) 

Q( fc ) 



s.t. tr (Q {k) ) < P {k) Vfc 

is A max (e^hW^ r pik ; 

Proof: Define R = Y!Li H^Hf fe) PW. Using Lem- 
mas Q] and |2l it is directly seen that A max (R) < A max (llj . 
showing the achievability of this solution, let 
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with v = v max (R). Apparently, tr (Q^) < P^ k) is satisfied 
and we can write 

A max (R) =v max (R) H ^H^qWhW^J v max (R) (5) 



K 



> 



(k) 



H 



(6) 



\k=l 



= J2 v H H( fc )H( fe ) ff p( fc )v = A max (R) , (7) 

k=l 

where © follows from the Rayleigh quotient and (0 is 
obtained by using d?} and writing v and v H inside the sum. 
Thus, choosing Q( fc ) as in (|4|i leads to the optimal value 

Amax ( R^ ■ H 

With this optimal choice of the covariance matrices the 
achievable sum-rate for joint relaying can be expressed as 



i? s =log 2 (l + ||h|| 2 P r 



log 2 



|h|l 2 P, 



1 + A max ( R 




(8) 



IV. TDMA-Based Relaying 



In this section, we will introduce a relaying scheme based on 
TDMA. This scheme includes a division of the transmission 
in K time slots, where user k occupies the k-th time slot 
exclusively. Also the relay incorporates this slot structure, i.e., 
the relaying matrix ¥^ k ' in time slot k can be adapted to 
the channel of user k only. Thus, the TDMA slot structure 
decomposes the channel in K independent single user relay 



A. Single-User Relaying 

In order to describe a single-user relay channel with a 
consistent notation, we assume the same channel model as 
introduced in subsection III-BI and assume there is only K = 1 
transmitting user. For this scenario, it was shown in [3] 
Theorem 1], that the optimal structure of the source covariance 
matrix Q^ 1 ' and the relay matrix are given by 

q(d = vWe^vW* 

F« = UftEfuW* 
= UWE«VW H is the SVD of H^, 



where 

£(*) 



diag(o^ 1} ,..., 



min{MW,M r } 



,1/2, S 



(1) 



diagfaW . . . , q%,\ and E, = diag(/f \ f^f/\ 

With these expressions, the achievable (sum) rate of 
the only user can be derived by evaluating ||2), which, together 
with the power constraints (Q]) allows the formulation of an 
optimization problem (cf. O): 



n 
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f (i) 

• > 3 M r 
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R^ 1 



N / 2 (i) (i) f (iy 
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offP + 1 



i=l 



(1) 



E ^ < p 

i=l 



qV>0 Vi, 



where N = min {N, M r , ikfW }. 

In the case of N = 1 receive antenna considered here, the 
problem is even simpler to solve. As in the joint relaying 
scheme, we have a\ — ||h|| 2 and <7 2 = . . . = a n = 0. 
Thus, the optimal choice of F is again obtained by choosing 
h = PrO- + a^q^)- 1 , h = ■■■ = fM r = 0. Also for the 
transmit covariance matrix, only q^ has a contribution to the 



rate, which makes q[^ = PW, g 2 = . . . = q M ' (1 ) = the 
optimal solution. Using this solution, the achievable rate i?W 
can be written as 



R(V = log 2 1 



||h|| 2 4 1) P( 1 )P r 



|h|l 2 R. 



a{ 1] PW 



log 2 (l + ||h|| 2 P r ) -log 2 1 



|h|| 2 P r 



1 + 4%', 



B. TDMA-based Transmission Scheme 

In order to decompose the if -user MARC in K single- 
user relay channels, we use a TDMA scheme, such that 
user k transmits only in a time slot of duration ' with 
X^fcLi T ^ = 1- I n eacn tmie s l ot ' me optimal choice of the 
relay matrix FW and the transmit covariance matrix Q' fe ) can 
be obtained as in subsection IIV-AI The only difference is the 
transmit power constraint: As user k only transmits in r^ k ' 
fraction of the time, it can use a transmit power of P' 1 '/t'" 
and still fulfills the average transmit power constraint (Q~|i. 
Thus, the rate of user k is given by 



R (k) 



log 



l+llh|MP. 



log 2 1- 



IhlPR. 



1 



(fc) ptk) 



and the sum-rate can be calculated as -Rsjdma = 
J\ = i-R^- This sum-rate can be optimized by the choice 
of tW, . . . ,t( k \ where the optimum is given as follows: 

Theorem 4: For the considered TDMA transmission 
scheme in the if-user MARC, the optimal durations of the 



time slots 



-{K) 



are given by 



r (fe) = 

opt 



K 



(9) 



pU) 



Proof: Consider the optimization problem 



max i? S ,TDMA(T") 

r(i),...,r< A ") 



s.t. h(r) = 1 - r(fc) =°. 

k=l 

where we introduced the vector r = [t^ , . . . , r^l and 
wrote i?s.TDMA as a function of r to be more conform with the 
notation used in optimization theory. First, it is easy to show 
that the above problem is convex and thus the Karush-Kuhn- 
Tucker (KKT) conditions provide necessary and sufficient con- 
ditions for optimality. For the underlying problem, the KKT 
conditions of a solution t* to be optimal can be formulated 
as h{r*) = and Vi? s ,TDMA("r*) + v*Vh(r*) = 0, where 
v* G K can be chosen arbitrarily. If we choose r* as suggested 
in Theorem [4] the first condition is obviously fulfilled. For the 
second condition we can calculate the derivations of -Rsjdma 
and h at the point t = t* as 



dR 



SJDMA 



<9r( fc ) 



= lo 



82 



1 



|h|l 2 R 



5(1 



log 2 (e)5||h|| 2 P r (l 



|h||2P r ). 

l|h 2 ||Pr) 



(S(l + ||h||2P r ) + l) 2 (l 



(S(l + ||hpP r ) + l) 



dh 



= -1, 



T—T* 

K 



respect to r' fc - ) are the same in all components k = 1, . . . , K, 



fe=i 



where S = ( ^ ) . As the derivations with 



we can select v* = 



a T (fe) 



and also the second KKT 



condition is fulfilled. 
Using this optimal time slot durations, the achievable sum-rate 
of the TDMA-based transmission scheme is given by 

/ \ 



Re 



TDMA 



= log 9 l + ||h|| 2 P r 



log 2 



1- 



V 



1+ J2 a^PM 



I 
(10) 



C. Comparison with Joint Relaying 

The comparison between the proposed TDMA-based and 
the previously considered joint relaying scheme is given by 
the following theorem: 

Theorem 5: In the considered K-user MARC, the TDMA- 
based transmission scheme with optimal time slot durations 
tW achieves at least the same sum-rate as the joint relaying 
scheme, i.e., 

-Rs,tdma > -Re- 

Proof: Comparing the rate expressions of -Rs.tdma and 
i?s from ( [Tol l and (0, respectively, it can be seen that it is 
sufficient to show that Y^f=i a i j) pU) > A max (ft) - This is 
easily obtained by writing 

f>Wp« = VA^fH^HWV^ (11) 



3=1 



3=1 



> An 



(12) 



where (TTZt follows from the convexity of the maximum 
Eigenvalue operation. ■ 
Note that (cf. |8 1) in (TTZb equality holds only if all Eigenvalues 

A max (ni j) ni j)H P( j n have the same Eigenvector. However, 

if the matrices Hf are statistically independent, which is 
a valid assumption in practical systems, this happens with 
probability zero. The consequence of this is that the TDMA- 
based approach is strictly better. 

A plausible explanation for this is the fact, that in joint 
relaying, the relay always has to come to a compromise, such 
that not only one but all incoming signals are amplified as 
good as possible. In the TDMA system, the relay has only 
one incoming signal, which can be optimally amplified by 
optimizing the relay matrix just for this signal. 

V. Simulation Results 

In order to confirm the result of Theorem[5] we will evaluate 
the achievable rates of the TDMA-based and joint relaying 
scheme by some simulations in this section. As stated in the 
preceding section, the TDMA-based scheme is strictly better 
if the channel matrices are statistically independent. 

However, the question for the quantitative gain of the TDMA- 
based scheme has to be investigated. 

For this purpose, we assume channels with independent 
Rayleigh-fading, i.e., all entries of the channel matrices 
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Fig. 2. Comparison of TDMA-based and joint relaying scheme for K = 8, 
MC 1 ) = . . . = MM = M r = 4, P« = . . . = P( 8 > = 10 

are ~ CJ\f(0, 1) and independent from each other. All results 
presented in this section are obtained by averaging over 1000 
channel realizations. 

For a system with K = 8 users, the achievable sum-rates are 
plotted in Figure [2] We assumed that all transmitters and the 
relay have M ( k > = M r = 4 antennas, while, as a prerequisite 
of this scenario, the receiver is assumed to have a single 
antenna. The transmit power of the users is assumed to be 
fixed at PW = 10 Vfc, while the power of the relay P r is 
varied. It can be observed that, especially for high values of 
P r , the performance of the TDMA-based scheme is up to 
10% better than with the joint relaying scheme. From further 
simulations, which are not visualized in this work due to the 
page constraint, it could be seen that this gain grows with the 
number of antennas at the transmitters M^ k \ the number of 
antennas at the relay M r , and with the number of users K. 

VI. Conclusion 

We have considered a K-user MARC, where all nodes 
except the receiver have multiple antennas. For this channel, 
we found the optimal solution for a previously introduced 
transmit strategy [7 1, where all transmitters access the channel 
jointly. Contrary to this approach, we also proposed a trans- 
mission strategy that is based on a TDMA structure, where 
each transmitter accesses the channel exclusively. As TDMA 
decomposes the if -user MARC in K single-user channels, 
the transmit strategies of single-user relay channels can be 
used. After optimizing these strategies and finding the optimal 
time slot durations for the scheme, it could be shown that the 
achievable sum-rates of the TDMA-based scheme are larger 
or equal than the rates achievable by the previous scheme. 
Moreover, equality of the rates only occurs for cases that have 
probability zero in practical systems. Finally, a quantitative 
sum-rate gain of around 10% was observed for one specific 
scenario by simulations. 
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Appendix A : Proof of LemmaQ] 

Let v max be the Eigenvector corresponding to the Eigenvalue 
Amax (B) with ||v max || = 1. Then, we have 

A m ax (B) = V^ ax Bv max 

< v ma X Bv max + V^ ax • (A - B) • V max 

V H A v 
v max -"- v max 

V max v max 

< A max (A), 

where the first inequality is due to the positive semidefiniteness 
of A — B and the last inequality can be derived from the 
Rayleigh quotient. ■ 

Appendix B: Proof of Lemma[2] 
Lemma [2] can be equivalently formulated as 

x H A(PI - U)A H x > Vxel™ 

Hence, it is sufficient to show that PI — II is a positive 
semidefinite matrix. Let Ai, . . . , A„ be the Eigenvalues of II 
with Ai > A2 > . . . > A„. Then, we have 

n 

P>tr(n)=^A. i >Ai, 

i=l 

where the second inequality holds because we have II >z and 
thus Ai > Vi. Using A m i n (PI — II) to denote the minimum 
Eigenvalue of PI — II, we can write 

A mm (Pi n) = p A max (n) = p - Ai > 0, 

i.e., PI II is positive semidefinite. ■ 
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